Heat capacity of suspended phonon cavities 
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We present a detailed analysis of the vibrational spectrum and heat capacity of suspended meso- 
scopic dielectric plates, for various thickness-to-side ratios at sub-Kelvin temperatures. The vi- 
brational modes of the suspended cavity are accurately obtained from the three-dimensional (3D) 
elastic equations in the small strain limit and their frequencies assigned to the cavity phonon modes. 
The calculations demonstrate that the heat capacity of realistic quasi-2D phonon cavities approach 
the linear dependence on T at sub-Kelvin temperatures. The behavior is more pronounced for the 
thinnest cavities, but takes place also for moderately thick structures, with thickness-to-side ratios 
7=0.1 to 0.2. It is also demonstrated that the heat capacity of the suspended phonon cavities is 
invariant under the product of the temperature (T) with a characteristic lateral dimension (L) of 
the sample. The present results establish a lower bound for the heat capacity of suspended meso- 
scopic structures and indicate the emergence of the quantum mechanical regime in the dynamics of 
bounded phonon cavities. 
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I. INTRODUCTION 

Suspended nanostructures have become relevant ele- 
ments for both basic research and technology. Micro- 
electromechanical systems (MEMS), such as cantilevers, 
gears and membranes, already find widespread use in 
several technological applications 1 . At the same time, 
current developments in surface nanomachining render 
possible the controlled fabrication of a large variety 
of suspended nanostructures 2 -! 3 - having, in particular, 
an extremely weak thermal coupling with the environ- 
ment. As a consequence, ultrasensitive bolometers 4 - and 
calorimeters^ have been developed with unprecedented 
sub-attojoule resolution, for operation in the T < 5 
Kelvin temperature range, envisaging the possibility of 
measuring the heat capacity of nano-objects and, even- 
tually, even single molecules^. In the realm of funda- 
mental research, the operation of nanoelectromechanical 
structures (NEMS) is finally approaching the quantum 
regime&SJii. The construction of suspended solid state 
quantum logic gate s 11 ' 12 is amid the applications antic- 
ipated for such structures, since the electron-phonon in- 
teraction, which is a source of decoherence and dissipa- 
tion for both quantum dot qubits^^ and single-electron 
transistors (SET>iS,, can be controlled in themi&±L±&±2. 

In fact, for most of the cases mentioned above the de- 
vices are operated at sub-Kelvin temperatures; the re- 
quirement of ultracold temperatures being specially se- 
vere for the operation of quantum logic gates. For in- 
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stance, in recent experimental realization, quantum dot 
charge-qubit a 13 : 14 and suspended SET— have been oper- 
ated at a base temperature of 20 mK. It is therefore rea- 
sonable to expect that suspended nanostructures com- 
prising such quantum devices as well as ultrasensitive 
bolometers and calorimeters will be functioning at tem- 
peratures T < IK. 

Despite the interest, there is not yet a comprehensive 
theory for the electron-phonon interaction in suspended 
nanostructure s 20 : 21 . A central issue of the problem is the 
difficulty of rigorously describing the low temperature 
acoustic phonon spectrum in suspended nano-devices, 
which is fundamental for determining (i) the electron- 
phonon interaction with its many consequences for the 
device operation and (ii) the device's thermal proper- 
ties, such as its thermal conductivity and heat capac- 
ity. At sub-Kelvin temperatures, the formalisms adopted 
for bulk materials may not produce correct results for 
the phonon spectrum of suspended nanostructures be- 
cause the wavelength and mean free path of the dominant 
phonons can be bigger than the physical dimensions of 
the structure. Moreover, there are no general analytical 
solutions for the vibrational modes of bounded suspended 
plates 2 ^. 

Motivated by these circumstances this work presents 
a detailed study of the phonon spectrum and the heat 
capacity (CV) of suspended rectangular dielectric nanos- 
tructures of various thicknesses at sub-Kelvin tempera- 
tures. The vibrational modes of the suspended cavity are 
accurately calculated from the three-dimensional (3D) 
elastic equations in the small strain limit and the ob- 
tained frequencies assigned to the cavity phonon modes. 
After obtaining a reliable phonon spectrum, with conver- 
gence assured for a few thousand cavity modes, the heat 
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capacity of isolated suspended mesoscopic phonon cav- 
ities having 3D and quasi-2D character is investigated. 
For such systems, the calculations demonstrate that the 
temperature dependence of Cy approaches the linear 
regime in sub-Kelvin temperatures, the effect being more 
pronounced for quasi-2D nanostructures. Nonetheless, a 
simple model of plane waves yields a phonon spectrum in 
good agreement with the 3D elastic model for the very 
thick suspended nanostructures. A dimensional analysis 
of the free vibrational modes also reveal that the heat 
capacity of the rectangular phonon cavities has the scale 
invariant form, that is, CV is invariant with respect to 
the product of the temperature with a characteristic lat- 
eral dimension. The present results indicate that the low 
temperature heat capacity of quasi-2D suspended nanos- 
tructures may have been underestimated and, therefore, 
sets a lower bound for their heat capacity. 



II. THEORETICAL FORMULATION 

We consider suspended nanostructures of rectangular 
geometry, with lateral dimensions defined by L x and L y , 
and thickness L z . Such a choice is motivated by the fact 
that several recent experiments 4 - ' 23 ' 24 , probing thermal 
and electrical properties of suspended nanostructures, 
have utilized square or rectangular plates with thickness- 
to-side ratios (7 = L z /L y ) including 0.1 [HI, 0.04 [H, 
and 0.015 The first structure is considered to be a 
moderately thick plate, while the last two cases are ex- 
amples of thin plates. Because the suspended nanostruc- 
tures are usually made of non-crystalline materials, like 
poly-silicon and amorphous SiN (Silicon Nitride), in this 
work the phonon cavities are taken to be homogeneous 
and isotropic rectangular structures. 

Because at sub-Kelvin temperatures the dominant 
phonon modes are long wavelength acoustic ones, with 
a mean free path that exceeds the dimensions of the 
structure^ ' 23 ' 24 , we resort to the elasticity theory to ob- 
tain the phonon spectrum of the suspended nanostruc- 
tures. In this limit the phonons correspond to the free 
vibrational modes of the cavity 1 ^, as determined by the 
elastic theory of solids 2 ^. The continuum elasticity model 
has been successfully used to describe the properties of 
propagating phonons in beam s 26 ' 27 , thin membrane s 16 ' 28 , 
and arrays of nanomechanical resonator s 29 ' 30 . 

To secure the correct description of the thermal prop- 
erties of suspended nanostructures of a few /im 2 in area 
in the sub-Kelvin temperature regime, at least a few 
thousand vibrational modes have to be calculated with 
confidence. A variety of methods intended to calculate 
the free vibrations of thick plates have been developed 22 . 
Since the simplified models like the Classical Plate The- 
ory (CPT)2i are adequate only for the lowest modes of 
thin plates, a three-dimensional analysis of the free vi- 
brations of the cavity is necessary. In general such meth- 
ods utilize the Rayleigh-Ritz formalism to determine the 
displacement field, which is represented as a series of or- 



thogonal polynomial s 32 ' 33 . In this work we follow the 
procedure developed by Zhou et alM- due to its simplic- 
ity and generality as well as for producing very accurate 
natural frequencies. For the sake of completness, the 
method is summarized next. 

For the problem of the free vibrations of isotropic 
structures in the small strain approximation, the kinetic 
(T) and strain (U) elastic energy functionals for the dis- 
placement field u(r, t) = XJ(r)e luJt can be written as 
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where p is the mass density, E is the Young's modu- 
lus and v the Poisson's ratio. The A quantities in the 
strain energy term are: Ai = J2i £ Hi A2 = J2i £ ii an d 
A3 = J2i<j £ iji f° r h3 — { x iUi z )i with the components 



of the strain E; 



diUi and £i 



djUi 



diUj. It is 



convenient to normalize the coordinates with respect to 
the dimensions of the plate, defining the dimensionless 
variables £ = 2x/L x , rj = 2y/L y and ^ = 2z/L z in the 
interval [-1,1]- The time-independent displacement field 
U(r) is then written as a sum of orthogonal Chebyshev 
polynomials multiplied by boundary functions Fs(l;,r]) 



l7x(£,»?,C) = F x (Z,r } )Y^A ijk P i (&P j (Ti)I\(C) (3) 

i,j,k 

Uy(tvX) = F y (z, v )J2BmPi(0PMPk(0 (4) 



i,j,k 



with the summations beginning from zero. The functions 
Pnix) are Chebyshev polynomials of the first kind and 
degree n, defined by the relation 



Pn(x) = cos [n arccos(x)] 



(6) 



with n a non-negative integer. The boundary func- 
tions F x (£, rf), Fy(£, rj) and F z (£, rj) have the general form 
F s(£,v) = fs(Ofsiv), with S = x,y,z. For our pur- 
poses the boundary conditions of interest are FF (free- 
free), CC (clamped-clamped) and CF/FC, which corre- 
spond, respectively, to the functions = fs(v) = 1; 

fHO = fi(v) = 1 - V 2 ^ and /KO = f}(v) = 1 ±V- 

Substituting the series representation for u(r, t) = 
U(r)e tuJt into Eqs. (P) and |2]), one obtains T max and 
U raal , which are the maximum values of T and U during 
a vibratory cycle. The frequency determinant is formu- 
lated by minimizing the functional U max — T max with 
respect to each of the coefficients {A}, {B} and {C}, to 
produce the 3D elastic equations of motion 
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TABLE I: Physical constants of relevant materials in the 
amorphous phase. 



Material 


p(g/cm J ) 


£(GPa) 


V 


GaAs — 


5.1 


71 


0.32 


Si^ 


2.3 


170 


0.22 


SiN^ 


3.1 


285 


0.20 


SiC-21 


3.0 


400 


0.20 



where Q = loL x ^J p/E is a dimensionless parameter and 
lo is the free vibration frequency, to be assigned to the 
phonons. In Eq. (O, [K] and [M] denote the symmet- 
ric stiffness matrix and the block diagonal mass matrix, 
respectively, which can be found in explicit form in Ref. 
[32| . Another useful dimensionless parameter associated 
with the frequency is A = (f2/A77r 2 ) y/l2(l — u 2 ), with 
A = L x /L y and 7 = L z /L yi which yields the normalized 
frequencies for the family of all rectangular plates with 
the same aspect ratios (A and 7) and elastic constants 
(E and v). 

An important aspect of the present analysis is the relia- 
bility of the phonon spectrum to be used in the evaluation 
of the thermodynamical properties of the nanostructure. 
For that reason the convergence of the highest frequen- 
cies was limited to be within 5%. The spectrum span can 
be extended by increasing the amount of basis functions 
P n (x) used in the representation of U(r). In the in- 
stance of the thickest square plate to be considered, with 
7 = 0.5, we have used n x — n y — 29 and n z = 13 Cheby- 
shev polynomials, yielding approximately 12000 reliable 
modes. For the thinnest plate, 7 = 0.02, the best results 
were obtained for n x — n y — 51 and n z — 4 that yielded 
4000 reliable phonon modes. The number of reliable fre- 
quency modes will determine the maximum temperature 
(T max ) for which the heat capacity can be calculated with 
confidence. 

The dependence of the phonon frequencies on the ma- 
terial parameters is such that higher frequencies are ob- 
tained for stiff and light materials. Table [T] contains the 
values of the mass density, Young's modulus and Pois- 
son's ratio for materials of relevance for the fabrication 
of NEMS. This work investigates nanostructures made 
of amorphous silicon carbide (a-SiC) because of its high 
rigidity and widespread use in the fabrication of sus- 
pended NEMS. 

The calculated spectra exhibit significant dependence 
on the dimensions of the nanostructure, as illustrated 
in Figure [T] for the first 2500 vibrational modes of free 
standing a-SiC mesoscopic structures. The structures 
have the same lateral dimensions L x — L y — L — 2 pm 
but different thickness-to-side ratios: 7= 0.02 (solid), 
0.05 (dashed) and 0.1 (dot-dashed). At the lower part 
of the spectrum the frequencies are higher for the thick 
plates, however, the behavior is reversed as the mode 
index a increases. The frequencies are also inversely 
proportional to the area of the plate. Moreover, by 
a numerical analysis of the eigenfrequencies it was 



observed that the vibrational spectrum of the cavities 
can be very well described by the form lo = loqoi^ in 
two limiting cases: for the quasi-2D phonon cavities 
(7 < 0.02) the fitting yields <f> < 1, whereas for the 
three-dimensional (thick) phonon cavities (7 > 0.2) one 
obtains 0.4 < <fi < 0.5 . Between the two cases, i.e. for 
moderately thick nanostructures, the frequencies cannot 
be well described by a single power curve. 




FIG. 1: Natural vibration frequencies, as a function of the 
mode index a, calculated by the 3D method for a free standing 
square (A = 1) a-SiC nanostructure of sides L = 2pm and 
7 = L z /L equal to 0.02 (solid), 0.05 (dashed) and 0.1 (dot- 
dashed) . 



In the following we examine the heat capacity of 
the phonon cavities, as predicted by the 3D analysis. 
For the sake of comparison, we also utilize a basic 
model to describe the confined phonons, which comprises 
some of the reductionist features commonly found in the 
literatur e) 7 ' 38 . Its main assumption is that the phonons 
can be described by plane waves with three independent 
polarizations: one longitudinal and two transverse. In 
addition, as the dimensions of the structure become suf- 
ficiently small, i.e. comparable to the mean free path of 
the phonons, these become standing waves satisfying the 
appropriate boundary conditions. The method is here 
designated Bounded Plane Wave Model (BPWM). For 
the nanostructures under consideration and because of 
the very low temperature, it is assumed that the phonons 
form standing waves in all three directions. In compar- 
ison with the 3D analysis, different predictions for the 
specific heat are expected on the basis of the BPWM, 
owing to its naive representation of the phonon modes in 
suspended nanostructures. However, despite the simplic- 
ity, it will be shown that the BPWM can describe the 
heat capacity of thick phonon cavities quite well. 

According to the BPWM, the phonon spectrum of a 
freely suspended nanostructure is easily obtained from 
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the wave vectors 



k 2 

^Imn 



(8) 



with I, m and n integers. The frequencies for the longitu- 
dinal and transverse modes are given by ui 



and 



u> 



t 

hnn 



I 

hnn 



V-Klr, 



V-Kl r , 



with the sound velocities of the a- 



SiC obtained from the elastic constants of the material: 
v l = 12, 170 m/s and v* = 7, 450 m/s. 



III. HEAT CAPACITY OF SUSPENDED 
PHONON CAVITIES 

Having calculated the displacement modes U Q (r) and 
the associated eigenfrequencies uj a , corresponding to the 
free vibrations of the plate, the quantum mechanical 
phonon modes of the cavity are obtained by the stan- 
dard quantization procedure^. As a result, we ascribe 
the energy U = *^2 a {n a + l/2)hw a to the phonon system 
and calculate the quantum mechanical heat capacity of 
the phonon cavity as 



C V (T) 



dU 
df 



d_ 

df 



E 



hllJn 



exp (hio a /kBT) — 1 



(9) 



with n a given by Planck's distribution. Since the har- 
monic regime complies with the small strain limit that is 
assumed in the present derivations, the constant-volume 
specific heat (c„) must equal the constant-pressure spe- 
cific heat (c p )22. Moreover, for bulk metallic samples it 
is generally found that the low temperature specific heat 
varies as c v — AT + BT 3 , comprising the electron and 
phonon contributions, respectively. 

Next we show predictions for the heat capacity ob- 
tained through the 3D analysis as well as results gained 
by the simplified method. The phonon cavities to be 
considered are free standing a-SiC square plates (A = 1) 
with lateral dimensions L = 2 /im and having different 
thickness-to-side ratios, namely 7 = 0.02, 0.05, 0.1, 0.2 
and 0.5. Cy is calculated for temperatures T ^ T max , 
where T max is the maximum temperature that allows re- 
liable results to be obtained with the available phonon 
modes. That is, if T > T max additional modes must be 
included in the calculation of Cy since the occupation of 
the high energy modes is increases. Figure [Ua) presents 
the specific heat c v — Cy fV as a function of temperature 
for phonon cavities of different thicknesses, as obtained 
through the 3D analysis. For temperatures T < 100 
mK the calculations reveal more than an order of mag- 
nitude difference between the quasi-2D (7 ^ 0.05) and 
the three-dimensional (7 > 0.2) phonon cavities. Figure 
[2Kb) shows the heat capacity (CV(T)) for two limiting 
cases, represented by the quasi-2D (7 = 0.02) and fully 
three-dimensional (7 = 0.5) suspended nanostructures. 
Results gained from the BPWM are also shown by the 
gray curves. Particularly, the BPWM predicts very good 
results for the three-dimensional cavities, but seriously 
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FIG. 2: a) Specific heat (c„ = Cv/V) as a function of temper- 
ature obtained from the 3D analysis for a free standing a-SiC 
square cavity of lateral dimensions L = 2 fim and 7 = 0.02, 
0.05, 0.1, 0.2 and 0.5, in that order from top to bottom, b) 
The heat capacity (Cv) for the quasi-2D 7 = 0.02 (solid) and 
fully 3D 7 = 0.5 (dashed). The gray curves are predictions 
from the BPWM (refer to the text). 



underestimates the heat capacity of the quasi-2D struc- 
tures. 

Throughout the analysis we have considered thin as 
well as thick suspended nanostructures. That raises the 
question, how is the system's dimensionality reflected on 
the behavior of Cv{T)l In Ref. [38| the dependence 
of the specific heat on the dimensionality of the system 
was investigated with a model similar to the BPWM. It 
was shown that the relation Cv oc T d should hold for a 
confined phonon gas in the low temperature limit, with 
d as the system's dimensionality, supporting the inaccu- 
rate notion that Cv oc T 2 for quasi-2D phonon cavities 
at sub-Kelvin temperatures. Here we perform such an 
analysis and demonstrate instead that the heat capacity 
of realistic quasi-2D phonon cavities approaches the lin- 
ear dependence Cv oc T in the low temperature limit. 
For that purpose consider the quantity 



Pc(T) = T 



d{\nCy{T)) 
dT 



(10) 



which provides the temperature dependence of Cy(T). 
For instance, if the heat capacity is given by Cy oc T Q , 
we have simply pc = ex. 

Figure[3]presents the calculated values of pc{T) for free 
standing square plates with lateral dimensions L = 2/xto 
and the thickness-to-side ratios previously considered, 
ranging through the quasi-2D to the fully 3D cases. The 
upper panel shows that Cy approaches the linear Cy oc T 
behavior at temperatures T < 0.3 K, particularly in the 
case of the thinnest cavities with 7 = 0.02 and 0.05. Such 
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0.1 0.2 
T(K) 



ral modes, which reduces the dimensions of the problem 
from three to two by incorporating some of the plate's 
characteristics such as bending moment a 25 ' 31 , yields a 
close estimate for the temperature dependence of Cy(T). 
The CPT fails, however, as the temperature raises above 
T > hirv s j '{L z ks) , v s being the sound velocity, because 
longitudinal and torsional modes begin to contribute sig- 
nificantly to Cy . Because of the significant difference be- 
tween the predictions made by the 3D analysis and the 
BPWM for thin nanostructures, pc may be a convenient 
observable to experimentally determine the emergence of 
coherent quantum mechanical dynamics in mesoscopic 
phonon cavities. 

An additional property of the suspended phonon cavi- 
ties is the scale invariant character of their heat capacity, 
described as Cy = ^(TL), where T represents the func- 
tional in Eq. Namely, Cy is invariant regarding 
the product of the temperature (T) with a characteristic 
lateral dimension (L) of the structure. For the sake of 
clarity we consider a square cavity, but the same result 
can be derived for rectangular cavities with L x — XL y , or 
triangular ones. First notice that the phonon frequency 
can be written as 



FIG. 3: pc(T) for free standing a-SiC square plates with L = 
2 //m. In the upper panel the predictions of the 3D analysis 
for different 7, as indicated by the labels. In the lower panel 
pc for the quasi-2D (7 = 0.02) and fully 3D (7 = 0.5) cases. 
The gray curves are the results gained from the BPWM. 



a result is expected to hold for strict 2D systems like 
graphene, although the effect has been predicted also 
for the specific heat of a 2D array of nanomechanical 
resonators 2 ^ In fact, the sub-T 2 behavior is observed 
here even in the case of the moderately thick cavities with 
7 = 0.1 and 0.2, at lower temperatures. As the tempera- 
ture increases pc tends to 3, indicating that the thickness 
of the cavity becomes much larger than the dominant 
phonon wavelength; on the other hand, for vanishing 
small temperatures, lower than the fundamental vibra- 
tional energies, the phonon cavity behaves as a 0D sys- 
tem. In this case the specific heat decreases exponentially 
and pc diverges with T _1 . The Cy oc T behavior is com- 
monly associated with quasi-lD systems like single-wall 
nanotubes (SWNT) at low temperatures^. However, as 
the temperature decreases beyond T < 0.5K a sublinear 
behavior is observed n such systems^. The effect is as- 
cribed to the overwhelming contribution of the flexural 
modes to Cy, since those modes present u)(q) oc q 2 and 
consequently Cy oc T 1 / 2 jHH. 

The lower panel of Figure [3] compares the results ob- 
tained from the 3D analysis with those of the simple 
BPWM. According to earlier calculations, both methods 
produce similar results for fully three-dimensional struc- 
tures, but the BPWM yields the wrong pc ~ 2 value for 
quasi-2D structures. In the last case, it has been verified 
that the simple classical plate theory (CPT) for flexu- 



7T 2 7 / E . , 

Ua = — A "V l2p(l-, 2 ) ■ 

Thus cu a oc 1/L, for plates of a given thickness-to-side ra- 
tio 7. The dimensionlcss parameter A a is also a function 
of 7 and A, therefore independent of the absolute dimen- 
sions of the plate. Then, from the definition of the heat 
capacity, Eq. (0), with the derivative and summation 
operations commuted, it is easily verified that a trans- 
formation that leaves oj a /T oc 1/(LT) invariant does not 
change the heat capacity. Consequently, the results that 
have been presented for square phonon cavities of lateral 
dimension L — 2/ito can be generalized for congruent cav- 
ities of arbitrary size L' , with the temperature re-scaled 
to T = (L/L')T. 

Different types of suspended cavities were also in- 
vestigated, such as bridges (CCFF) and cantilever-like 
(CFFF) structures, yielding results in qualitative agree- 
ment with those previously illustrated, for both the vi- 
brational spectrum and the specific heat. It is observed 
that the parameter pc shows a tendency towards the 
value 1 for elongated structures. For instance, in the 
case of 7 = 0.1 and A = L x /L y =4we obtained p c ~ 1.5 

at T — ^max* 

The hitherto calculations of the heat capacity of 
phonon cavities have not taken into account the addi- 
tional degrees of freedom comprised by impurities, disor- 
der and surface defects, etc. that will be responsible for 
an increase of Cy. In particular, the specific heat of bulk 
noncrystalline solids exhibits an anomalous linear varia- 
tion with the temperature^ for T < 1 K. The present 
results, however, set a lower bound for the specific heat 
of such dielectric nanostructures. 
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IV. CONCLUSIONS 

We presented a detailed investigation of the vibrational 
spectrum and the heat capacity of suspended dielec- 
tric mesoscopic structures of various thicknesses at sub- 
Kelvin temperatures. More than 4000 frequency modes 
of the cavity were accurately obtained from the 3D elas- 
tic equations in the small strain regime. It is therefore 
demonstrated that the low temperature heat capacity of 
realistic quasi-2D phonon cavities have an approximate 
linear dependence on T, a result that contradicts esti- 
mates obtained by simple models. The sub-T 2 variation 
of the heat capacity is observed even for the moderately 
thick mesoscopic structures. The results show the im- 
portance of a fully 3D analysis based on the elastic equa- 



tions of suspended plates, bridges and cantilevers, for the 
correct determination of their thermal properties. Fi- 
nally, the sub-T 2 effect evidences the quantum mechan- 
ical nature of the phonon cavity dynamics and sets a 
lower bound for their specific heat. The reported results 
should have special interest for suspended nanostructures 
intended to be part of solid state quantum devices. 
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